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Abstract 

The Brans-Dicke (BD) theory of gravity is taken as a possible theory of k-essence coupled 
to gravity. It then has been realized that the BD scalar field does indeed play a role of a k- 
essence, but in a very unique way which distinguishes it from other k-essence fields studied in the 
literature. That is, first in the BD scalar field-dominated era when the contribution from this 
k-essence overwhelms those from other types of matter, the BD theory predicts the emergence of a 
yet-unknown zero acceleration epoch which is an intermediate stage acting as a "crossing bridge" 
between the decelerating matter-dominated era and the accelerating phase. Upon realizing this, 
next, closer study of the effects of this k-essence on the evolutionary behavior of the matter- 
dominated and the accelerating eras has been performed. The result of the study indicates that 
the BD scalar field appears to interpolate smoothly between these two late-time stages by speeding 
up the expansion rate of the matter-dominated era somewhat while slowing down that of the 
accelerating phase to some degree. Thus with the newly found BD scalar field-dominated era in 
between these two, the late-time of the universe evolution appears to be mixed sequence of the 
three stages. 



* e-mail : hongsu@astro.snu.ac.kr 
t present address 



1 



I. Introduction 



Perhaps one of the greatest challenges in the theoretical cosmology today would be to 
understand the emergence and nature of the observed late-time acceleration of the universe 
and provide an answer to the so-called "cosmic coincidence conundrum" which concerns the 
puzzle : why at the present epoch are the energy densities of dark energy and of dust-like 
dark matter of the same order of magnitude ? For instance, according to the recent WMAP 
data |l| , the universe energy density appears to consist of approximately 4-per cent of that of 
visible matter, 21-per cent of that of dark matter and 75-per cent of that of dark energy. Up 
until now, the most conservative candidate for the dark energy is the cosmological constant A 
and perhaps the most fashionable candidates with non-trivial dynamics are quintessence and 
k-essence. And the main difference between the two models is that the standard quintessence 
models 2 involve canonical kinetic terms and the sound speed of c 2 s = 1 while the k-essence 
models 3[ employ rather exotic scalar fields with non-canonical (non-linear) kinetic terms 
which typically lead to the negative pressure. And the most remarkable property of these k- 
essence models is that the typical k-essence field can overtake the matter energy density and 
induce cosmic acceleration only at the onset of the matter-dominated era for some period, 
at about the present epoch. Usually the emergence of such exotic type of scalar fields with 
non-canonical kinetic terms has been attributed to the string/supergravity theories in which 
the non-linear kinetic terms generically appear in the effective action describing moduli and 
massless degrees of freedom. We also note that a pioneering work on accelerating cosmology 
which has some relevance to the current intensive study of scalar field models for dark 
energy had been presented by Peebles and Ratra some time ago. In this regard, here we 
focus particularly on the Brans-Dicke (BD) theory of gravity [5] (which is one of the simplest 
extensions/modifications of Einstein's general relativity) as it involves probably the simplest 
form of such non-linear kinetic term for the BD scalar field. Besides, the BD scalar field 
(and the BD theory itself) is not of quantum origin. Rather it is classical in nature and 
hence can be expected to serve as a very relevant candidate to play some role in the late-time 
evolution of the universe such as the present epoch. Thus in this work, we would like to take 
the BD theory of gravity as a possible theory of k-essence coupled to gravity and explore the 
role played by the BD scalar field in relation to the unified model for dark matter and dark 
energy. Therefore in our present study, the BD theory is viewed as a theory of k-essence 
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field (with non-canonical kinetic term) coupled non-minimally to gravity. Namely in this 
context, unlike the "scalar-tensor theory" spirit of the original BD gravity, the BD scalar 
field is not viewed as a (scalar) part of the gravitational degrees of freedom but instead is 
thought of as playing the role of a k-essence which is a matter degree of freedom. 

II. Brans-Dicke scalar field as a unique k-essence 

The spirit of Brans-Dicke extension of general relativity is an attempt to properly incor- 
porate both the Mach's principle Q and the Dirac's large number hypothesis (| in which 
Newton's constant is allowed to vary with space and time. In general, the Brans-Dicke 
theory of gravity, in the presencs of matter with Lagrangian Cm, is described by the action 

s = I ^ [ii - + C A (1) 

where $ is the BD scalar field representing the inverse of Newton's constant which is allowed 
to vary with space and time and u is the generic dimensionless parameter of the theory. 
Extremizing this action then with respect to the metric g^ v and the BD scalar field $ yields 
the classical field equations given respectively by 

G^v = Rpv ~ -nQnuR = 8ttT^ u d + -q-T™, (2) 



where 

rpBD 
flV 



1 

8^ 



^(V^V,$ - i^V Q $V Q $) + ^(V M V,$ - ^V Q V a $) 



T%, = p M9^ + (Pm + Pm)UJJ v 



and 



V « V "*-(^3) T " A -(2^3)(^'-^- < 3 > 

Note that = (2/ ' ^fg)b{^fgH M )lbg ll,v ', but T^ U D cannot be defined in a similar manner. 
This is due to the presence of the non-minimal coupling term ~ ^fg^R in the action which 
obscures the separation of the scalar ($) Lagrangian from the tensor {g^ v ) Lagrangian. 
Indeed, the expression for Tz, and the coefficient factor (87r/$) in front of originate 
from this non-minimal coupling term ~ ^fg^R upon extremizing the total action S in 
eq.(l) with respect to g^ u . Next, Pm and pm are the pressure and the energy density of 
matter consisting of those of radiation (denoted by rad) and dust (denoted by m), i.e., 
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Pm = Prad + P m = Prad and pu = Prad + Pm respectively. In this perfect fluid approximation 
for matter, U a = dX a /dr (where r denotes the proper time) is defined to be the 4- velocity 
of a perfect fluid element normalized such that U a U a = — 1. Note that the Einstein's general 
relativity is the u — > oo limit of this BD theory. Note also that in the action and hence 
in the classical field equations, there are no direct interactions between the BD scalar field 
$ and the ordinary matter. Indeed this is the essential feature of the BD scalar field $ 
that distinguishes it from "dilaton" fields in other scalar-tensor theories such as Kaluza- 
Klein theories or low-energy effective string theories where the dilaton-matter couplings 
generically occur as a result of dimensional reduction. (Here we would like to stress that we 
shall work in the context of original BD theory format not some conformal transformation 
of it. That is, we shall work in the Jordan frame, not in the Einstein frame.) As a matter 
of fact, it is the original spirit |5| of BD theory of gravity in which the BD scalar field $ is 
prescribed to remain strictly massless by forbidding its direct interaction with matter fields. 
In what follows, however, we shall first maintain the BD scalar field alone and drop all the 
others (radiation + matter) in order to uncover the role of the BD scalar field particularly in 
connection with the dark matter - dark energy mixture at the present stage of the universe 
evolution. In this context, the theory can now be viewed as that of a exotic scalar matter 
field coupled to gravity but not a pure BD gravity theory in which both the scalar (the BD 
scalar field $) and the tensor (the metric g^) fields belong to the gravitational degrees of 
freedom. Then one consequence of this is the fact that the statement such as 11 uj — > oo 
amounts to the Einstein gravity limit" loses meaning and the BD parameter uo remains a 
(yet) completely undetermined parameter of the theory. Later, however, we shall see that 
this BD parameter uj will be determined once we impose the energy-momentum conservation 
V^T^ = which is nothing but the Bianchi identity, i.e., a consistency condition. 
We now work in the spatially-flat (k = 0) Friedman- Robert son- Walker (FRW) metric by 
assuming the homogeneity and isotropy of the spacetime, i.e., 

ds 2 = -dt 2 + a 2 (t)\ , - +r 2 (d9 2 + sin 2 9d4>% 
1 — kr z 

$ = $(t) (4) 

where k = for the spatially-flat case and a(t) denotes the scale factor. Then in terms 
of this FRW metric, the Friedmann equation representing the Einstein equation and the 
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Euler-Lagrange's equation of motion for the BD scalar field read, respectively 



a J 



uj ( <E> 
6 V $ 



$ + 3 ( - ) $ = 0. 

.a, 



(5) 
(6) 



Thus in order to solve these coupled non-linear field equations, we start by assuming that 
the solution ansatz is given by 



$(t) = —a n (t), a{t) = a(0) 
Go 



1 + 



P 



/Go 



(7) 



where Go denotes the value of Newton's constant at t — from which the effective New- 
ton's constant, i.e., the inverse of the BD scalar field deviates and f3 is some dimensionless 
parameter. Note that (neglecting all matter contents) Go is the only fundamental scale in 
this theory. First, substituting the first of this solution ansatz into the Friedmann equation 
in (5) yields the algebraic equation for the power index n as um 2 — 6n — 6 = which in turn 
gives n = (3± + 9)/cu. Next, substituting the same ansatz this time into the BD scalar 
field equation in (6) leads to a non-linear equation for only the scale factor 

Finally, substituting the second of the solution ansatz into this equation allows us to deter- 
mine the power a as a = l/(n + 3). Therefore, a solution (which indeed turns out to be a 
general solution) to the coupled Friedmann-BD scalar field equation reads 



a(t) = a(0) 



1 + 



P 



n+3 



(8) 



m = ^ra n (t) = i-a"(0) 



1 + 



P 



/Go 



n+3 



with n — (3± \/Qu) + 9)/uj as obtained above. Note, here, that if we demand the expanding 
universe condition, the value of uj and hence that of n should be such that (n + 3) > 0. As 
we shall see later on, this condition will indeed be met in the final form of the solution. 
Note that in addition, the field equations in eqs.(5) and (6) should be supplemented by 
the energy- momentum conservation, V V T^ D = 0. Generally speaking, if we ignore the 
inhomogeneities arising from the (linear) field perturbations, the BD scalar field can be 
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treated as a perfect fluid with stress tensor T^ U D = T^f /G = PbdQ^u + (pbd + Pbd)U^U u . 
Then, V^T^, = is given, in terms of the FRW-metric, by 



Pbd = -3 (^j (pbd + Pbd)- 



(9) 



Note also that this energy-momentum conservation is indeed a consistency condition origi- 
nating from the (geometric) Bianchi identity V^G^ = with G^ being the Einstein tensor 
given above in eq.(2). Therefore, we now need to check if the general solution in eq.(8) 
satisfies the energy conservation equation given in eq.(9), which is indeed a consistency 
equation. As we shall see in a moment, this energy conservation equation allows us to de- 
termine the value of the BD cu-parameter to a particular one. Thus to this end, we begin 
with the evaluation of the energy density and the pressure of the BD scalar field. Thus from 
rpBD (i _ p BD fiv _|_ (p BD + Pbd)U^U u , we can read off the energy density and the pressure as 



„ rpt rp p rpr rpQ rp<f> 

Pbd — — ±t — l tt-, *bd — J r — 1 e — 1 <h- 



(10) 



Note also that the pressure can be calculated using an alternative expression, Pbd = (T* + 
T tt )/3. Thus using eq.(10) and the expression for the energy-momentum tensor for the BD 
scalar field given in eq.(2), we get 



1 



Pbd = 



Pbd — 



16%G 
1 

167rG n 



£ 
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A\ a - 



£ 



(ii) 



^ +2 -1 



Next, from the expressions for exact solutions eq.(8), it follows that 

-(n+3) 

f I [ I I 
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'£ 



'a 
n \ - 



n 



n + 3 



a(t) 



n 



+ <„-i)(H) 



a(0) 

2" 



(12) 



3n 



{n + 3} 



a(t) 
a(0) 



-2(n+3) 



where (3 = /3/y/Go. Therefore from eqs.(9), (11) and (12), the energy-momentum conserva- 
tion equation becomes 



uon 



6 



87rG V n 
which, in turn, gives 



£ 



£ 



UO = — - 



8nG \n 



n=-2 



£ 



£ 



+ 



uon 



n 



£ 



(13) 



and hence finally the solution reads 



a(t) = a(0) 



1 + 



Gn 



(-TQ 



1 + 







1 -2 



(14) 



Note also that this set of values in eq.(13) is consistent with n = (3 ± \/Guo + 9)/oo that has 
been determined from the field equations. As mentioned earlier, the field equations leave 
the BD cu-parameter and the power index n undetermined but it is the energy-momentum 
conservation (which is a consistency condition) that finally determines their values. Having 
obtained the explicit expressions for energy density and pressure of the BD scalar field (being 
treated as a perfect fluid), we now attempt to determine the equation of state and the speed 
of sound. Firstly, using eqs.(ll), (12) and (13), we have 



Pbd 



16nG 



Q{3 2 



'ajty 
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-2 



P, 



BD 



16ttG 
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which yields the equation of state to be wbd — Pbd/ Pbd = — 2/3 2 /6/3 2 
for later use, note that (pbd + 3Pbd) = 0), namely, 

1 



P; 



BD 



~3 PBD 



(15) 

•1/3 (and hence, 
(16) 



meaning that the BD scalar field turns out to be a barotropic (i.e., constant w) perfect fluid. 
Next, the computation of the sound speed needs more careful treatment. Namely, using 
again eqs.(ll) and (12), it follows that 
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8ttG 
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2/3 3 



(17) 



'a(t)' 
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Therefore, the speed of sound in this BD scalar field fluid turns out to be 

2 _ dP BD _ dP BD /dt 2p 3 1 

~ " ~3 



(18) 



dp B D dpBD/dt -Q/3 3 
where we followed the spirit of definition for the effective sound speed of perturbations sug- 
gested by Ratra Namely, it is negative. Indeed, the negative definite sound speed squared 
generically signals the instabilities. Types of ordinary matter such as radiation or dust ex- 
hibit non-negative pressure and hence non-negative sound speed squared. But for rather an 
exotic type of matter like the k-essence, it is not so surprising to have negative sound speed 



squared (as long as it is a barotropic fluid) as it may well result from the negative pressure. 
Indeed, c 2 in our BD scalar field fluid is essentially the same as the equation of state wbd 
and hence is negative. We shall now be more specific in relating the negative definite sound 
speed squared to the instabilities of perturbations. Indeed, this issue of instability plays a 
key role in the so-called XCDM model for dark energy |9( as well. And to this end, we shall 
follow the general formulation given by Peebles and Ratra j^] and apply it to the present 
model of BD scalar field fluid. In the context of linear perturbation theory in which the den- 
sity perturbation is described by psD{t,x) =< Pbd (t) > +8pBD{t,%) (where < ... > denotes 
the background (mean) value), the equation of energy- momentum conservation, V 1 v T^d = 
yields Spbd = c 2 V 2 5pbd where c 2 = dPso/dpsD as usual. 

(1) If c 2 s > 0, this perturbation equation becomes a wave equation whose solution would be 
given, say, by 5pbd = 5 poe~ l ^ t+lk ' x and the usual positive definite sound speed squared (or 
the real value of sound speed) indicates the propagating mode of perturbation. 

(2) But if c 2 < 0, the solution to this perturbation equation would be given by 5pbd = 
5 poe ± ^ t+ ' lk ' x and the negative definite sound speed squared (or the imaginary value of sound 
speed) indiates the (exponentially) growing or decaying mode signaling the instability of per- 
turbation. Besides, since c 2 = dP BD /dp BD = — \u\ 2 /\k\ 2 < 0, this "imaginary" sound speed 
indicates that the increase (decrease) in density leads to lowering (growth) of pressure sup- 
porting the emergence of instability. 

Some ways to avoid the negative definite sound speed squared and hence the unstable inho- 
mogeneities when encountered have been discussed in the literature [a]. Among them, one 
very convincing argument is due to Steinhardt [lflj . It suggests that although the imaginary 
sound speed might mean accelerated collapse of inhomogeneities, such instability could be 
avoided at least at the subhorizon scale by taking into account the dependence of the sound 
speed on the wavelength characterizing the instabilities. 
Thus far, we have determined the behaviors 

a(t) ~ t, pbd 5777, Pbd = -\pbd- (19) 

a z {t) 6 

Based on this observation, we now discuss the role played by the existence of this scalar field 
dominated era in relation to the so-called "coincidence problem" asking why the transition 
from the matter-dominated era to the (observed) cosmic acceleration phase occurs only at 
the present stage. First, regardless of the type of matter (i.e., whether it is the ordinary 
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ones like radiation and dust or the k-essence such as the BD scalar field), one can easily 
realize that the strong energy condition can be represented by 

= 47rG [3P + p] > (20) 

where £ M denotes the future-directed timelike unit normal vector and we used p = T lJLV ^^ l/ 
and T£ = (3P — p). Next, in the standard cosmology in which ordinary type of matter 
(radiation or dust) is coupled to Einstein gravity, the time-time component of the Einstein 
field equation gives 

®=-!*V + ,) (2!) 

whereas in the present case of k-essence model in which the BD scalar field plays the role of 
the k-essence, the time-time component of the metric field equation reads 



a\ uj ( <E>\ 1/$ 



mj 3 v$y 3 v$ 

Now from eq.(ll) above, we have 

1 



(3P BD + Pbd) = 
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4 - U +6 UU +6 U. 



(22) 



(23) 



16ttGo 

Then using eq.(12) and by substituting the set of values uj = —3/2, n = — 2 determined 
above, the eqs.(22) and (23) yield 

W = g- £ ( 3P ^ + ^)- (24) 

Thus to summarize, for both the ordinary matter case and the k-essence case, the strong 
energy condition is always represented by the factor (3P + p) and it is its sign which al- 
ways determines whether the universe decelerates or accelerates. With this preparation, we 
now recall the physical quantities characterizing each era in the universe evolution which 
are summarized in TABLE I. The four rows represent radiation-dominated, matter(dust)- 
dominated, (BD) scalar field-dominated and the accelerating eras, respectively. Of course, 
the total energy density p and the pressure P consist of all the contributions coming from 
each component, radiation, dust, BD scalar field and some unknown entity (denoted by the 
subscript "acc" ) leading to the acceleration, namely 

P = Prad + Pm + PBD + Pace, (25) 
P — Prad + Pm + PbD + Pace 

9 



Prad ^ a 4 (t) ' ^Tdd 3 Prad 


a(t) -t 1 / 2 , (3P rad + p rad ) > 


deceleration 


P™ ~ a 3 (t)> Pm - 


a(f) ~ t 2 / 3 , (3P m + p m ) > 


deceleration 




a(t) ~t, (3Pbd + Pbd) = 


zero acceleration 


Pace = consi., P acc = -wpacc (-1 < u; < -1/3) 


a(t) ~ t m (m > 1), (3P acc + Pacc ) < 


acceleration 



TABLE I: Summary of cosmic evolution in the presence of the scalar field-dominated era. 



with P m = as it is approximated as a dust. The presence of the last (i.e., the most 
recent) era, namely the acceleration phase has been introduced based on the cosmological 
observation of the present large scale structure such as the anisotropy in CMBR and the 
luminosity of type la supernovae at high redshift which all suggest that the universe is 
currently undergoing cosmic acceleration and is dominated by dark energy component with 
negative pressure And we placed the (BD) scalar field-dominated era in between the 
matter-dominated and the accelerating eras based on the way pbd scales with the scale 
factor a(t), namely pbd ~ l/a 2 (t). Thus from the scaling behavior of the energy density of 
each epoch, it is straightforward to see that as the universe expands (i.e., as a(t) grows), 
p ra d rises and falls first, p m next, then does pbd and lastly p acc begins to dominate. Of 
course, the behavior of other physical quantities during the (BD) scalar field-dominated 
epoch also strongly suggests that it should be inserted between the matter-dominated and 
the accelerating eras. Thus with this new era being inserted, it is interesting to realize that 
now we are witnessing a grand picture of universe evolution in which, as the universe expands, 
or as time goes on, the energy density p dilutes more and more slowly, the pressure P keeps 
decreasing from positive value eventually toward negative one, the scale factor grows more 
and more rapidly and lastly, the strong energy condition (i.e., the sign of (3P + p)) moves 
from "yes" towards "no" leading to the transit of cosmic evolution from deceleration towards 
acceleration past zero acceleration during the newly-inserted (BD) scalar field-dominated 
era. Among others, therefore, it appears that the newly found presence of the (BD) scalar 
field-dominated epoch provides a picture of smooth transition from the decelerating matter- 
dominated era to the epoch of current acceleration as the acceleration is zero during this 
era. Moreover, the BD scalar field itself is known to possess a generic dark nature as it is not 
allowed to have direct interactions (couplings) with ordinary matter (radiation + dust) from 
the outset, namely at the Largangian level since otherwise it would violate the cherished 
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equivalence principle as has been originally pointed out by Brans and Dicke themselves. 
Indeed, one of the features of the present study (in which the BD scalar field is taken to 
play the role of a k-essence field) that distinguishes it from other k-essence models can 
be summarized as follows. It provides no direct mechanism for the arrival (or emergence) 
of the present cosmic acceleration phase. Instead, it presents a good reason or evidence 
explaining why the matter-dominated era is to be followed by an accelerating phase as 
the acceleration is "zero" during the BD scalar field-dominated era in between the two. 
Namely, it provides a crossing bridge between the matter-dominated and the accelerating 
eras that has been missing from the scene thus far. From the theoretical perspective, this is 
certainly a physically natural and meaningful prediction of our model for dark energy. On 
the observational side, however, there has been no evidence for such zero acceleration epoch 
to date. 

III. Effects of BD scalar field (a k-essence) on the late-time universe evolution 

Thus far we have ignored the possible contributions from other types of matter and 
concentrated on the role played by the BD scalar field, i.e., the k-essence in order particularly 
to explore the (BD) scalar field-dominated era when the contribution from the k-essence 
overwhelms all others. And as a result, we found that such (BD) scalar field-dominated era 
is a yet-unknown zero acceleration epoch that should be inserted in between the decelerating 
matter-dominated and the accelerating eras acting as a "crossing bridge" between the two. 
Upon realizing this, then, our next mission should be a closer study of the effects of the 
k-essence (i.e., the BD scalar field) on the evolutionary behavior of the matter-dominated 
and accelerating eras. Thus we now should include all the three components, the dust-like 
matter, BD scalar field and the entity that drives the current cosmic acceleration. For the 
sake of definiteness of our demonstration, this entity responsible for the cosmic acceleration 
shall be taken as the cosmological constant A. Therefore, this study of the effects of BD 
scalar field on the nature of matter-dominated and accelerating eras can be thought of as 
our attempt to build a "unified model" for dark matter-dark energy that is currently under 
intensive exploration in the theoretical cosmology. 

We now reconsider the action for the Brans-Dicke theory of gravity, this time in the presence 
of the cosmological constant A with mass dimension 4 as well as the dust-like matter (denoted 
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by m), 
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A + Cr 



(26) 



16tt V $ 

Again extremizing this action with respect to the metric g^ v and the BD scalar field $ yields 
the classical field equations given respectively by 



where 

nBD 
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[(3P m - p m ) - 4A] 



(28) 



With Pry 



(2w + 3) 

in the dust approximation of matter. As usual we work in the spatially-flat 
FRW metric given in eq.(4) in terms of which the Friedmann equation representing the 
Einstein equation and the Euler-Lagrange's equation of motion for the BD scalar field are 
given, respectively 

'd\ 2 8tt . , . uj /$\ 2 /d\ 



$ + 3 



57T 



\{p m - 3P m ) + 4A] 



(29) 
(30) 



m/ (2a; + 3) 

Once again, we stress that these field equations have to be supplemented by the consistency 
conditions. Namely, in addition to these classical field equations for the metric and the 
BD scalar field $, there is one more set of equations which are consistency conditions as 
they result from the geometric Bianchi identity V l/ G At!/ = 0. Thus from eq.(27) we have 



= V V {R 



[IV 



-sTR) = v„ 



$ y 



(31) 



Note first that, according to the original spirit of Brans and Dicke |5[, in order not to 
interfere with the successes of the equivalence principle, the BD scalar field $ is assumed 
not to enter into the equations of motion of ordinary (dust) matter so that obeys the 
usual conservation law V^X^ = which, in terms of the FRW metric, takes the familiar 
form 



Pm = -3 ( - ) (Prn + Prr 



(32) 
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where again P m = in dust approximation. Therefore, we are left with 

= SttA^V" - STr^pTT + ^G Q V v [P BD g^ + (p BD + P BD )UW] (33) 
which, again in terms of the FRW metric, becomes 

Pbd + 3 (£) {pbd + Pbd) = ^(p™ + A) (^j . (34) 

Thus in order to determine the natute of late-time universe evolution, one in principle 
has to attempt to solve the coupled, non-linear field equations given in eqs.(29) and (30) 
subject to the consistency conditions in eqs.(32) and (34). However, one cannot solve them 
with both the dust-like matter and the cosmological constant present as they carry distinct 
equations of state. In practice, therefore, we consider the two stages ; first, the matter- 
to-scalar field-dominated era (MATTER-TO-SCALAR) transition period and second, the 
scalar field-dominated-to-accelerating phase (SCALAR-TO-ACCELERATION) transition 
period. 

1. MATTER-TO-SCALAR transition period 

At this stage, the energy density of universe is assumed to be dominated by those of dust-like 
matter, p m and BD scalar field, p BB >. First of all, in order to solve the coupled Friedmann 
equation and the BD scalar field equation in eqs.(29) and (30) (with A absent for the case at 
hand), the dependence of the matter energy density p m on the scale factor has to be given. 
As usual, this can be determined using the equation of state for the matter-dominated era, 
P m = and the energy-momentum conservation equation (32) which is one of the consistency 
equations. The result is the familiar dependence, p m ~ l/a 3 (t), or more precisely, 

(35) 



r<2 



a(t) 



_a(0)_ 

Then the coupled field equations for this transition period becomes 
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+ T i - ~ i , (36) 



■■ „/0\ • 8lTK 1 , . 

+ W (5J+3)gR (37) 



where k = a 3 (0) and of course these field equations have to be supplemented by the remaining 
consistency equation 

K 1 / $ 
r 



p BD + i(l)(p BD + P BD ) = — (38) 
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We now start with the solution ansatz 

1 



$0) = TT t 1 + X*F , a(t) = a(0) [1 + *t] £ 



(39) 



Substituting these solution ansatz into the coupled field equations (36) and (37) above yields, 



2(u + l) 
a = — ;r, 7 



(3a; + 4) ' (3cj + 4) 
Thus basically the solutions behave as 



X 



4vr (3u + 4) 
G~ (2uj + 3) ' 



$(t) ~ t 2 /( 3 -+ 4 ), a (t) ~ t 2 (-+D/(3-+4) 



(40) 



(41) 



and they have actually been known p] for some time in the Brans-Dicke cosmology Note that 
in the standard Einstein gravity limit u — > oo in which the dynamics of the BD scalar field 
is washed out, i.e., $ — > 1/Gq, one recovers the standard behavior for the matter-dominated 
era, a(t) ~ t 2 / 3 . With these exact solutions at hand, next we consider the behavior of the 
BD scalar field (being treated as a perfect fluid) in this transition period. Recall, first that 
its energy density and pressure are given by eq.(ll) earlier, which have been obtained using 
eqs.(2) and (10). Thus by plugging the exact solutions for the case at hand, eqs.(39) and 
(40) into eq.(ll), we end up with 



Pbd 



r<2 



(5w + 6) 



a(t) 



-2/a 



Pbd - 7^2 

Lr 



(2cj + 3)(3w + 4) \a(0) J 
2(w + l) fa(t)Y 2/a 



(42) 



(2w + 3)(3cj + 4) \a(0) 

where a is as given in eq.(40). Here a causion needs to be exercised. Namely, these expres- 
sions for p BD and Pbd do not necessarily mean that p BD < and Pbd > 0. Indeed, it is 
the other way around if we demand the positive-definite energy density for the BD scalar 
field, i.e., pbd > 0. Then it immediately follows that Pbd < 0. To see this is indeed the 
case, note that demanding pbd > leads to uj < —3/2 or —4/3 < to < —6/5 which, in 
turn, indicates {2u + 3) (3a; + 4) > and {uj + 1) < and also (5u + 6) < 0. Therefore, 
clearly Pbd < 0. Besides, this condition pbd > (particularly uj < —3/2) also guarantees 
the condition of cosmic expansion a = 2{uj + 1)/ (3a; + 4) > 0. Next note that as a result, 
the equation of state of the BD scalar field perfect fluid in this transition period is given by 

P BD 2(uj + 1) 



W B D 



Pbd 



(5a; + 6) 



< 0. 



(43) 



14 



Namely, it is negative implying that already in this matter-to-scalar field-dominated era 
transition period, the BD scalar field behaves as a negative pressure component. In the 
earlier section in which we maintained the BD scalar field alone and ignored the contributions 
from other types of matter, it has been realized that the equation of state of the BD scalar 
field is wbd — Pbd/ Pbd = — 1/3 implying the emergence of zero acceleration epoch driven 
by the negative pressure. This negative pressure nature of the BD scalar field as a k-essence, 
therefore, appears to begin already in this transition period. Next, the total pressure in this 
transition period is given by P tot = P m + Pbd = Pbd < 0. Thus it is the BD scalar field 
that renders the total pressure at this stage negative and this, in turn, implies that the BD 
scalar field would therefore act to speed up the expansion rate of the matter-dominated era. 
This actually will turn out to be the case since a > 2/3 as will be shown shortly in eq.(51) 
later on. Lastly, note that these solutions in eqs.(39) and (40) indeed satisfy the remaining 
consistency condition in eq.(38). It can be checked in a straightforward manner using the 
expressions for p BD and Pbd given in eq.(42) and the exact solutions in eqs.(39) and (40). 
This consistency condition, however, does not determine the BD u;-parameter to a particular 
value. 

2. SCALAR-TO- ACCELERATION transition period 

At this stage, the universe energy density is assumed to be dominated by those of BD 
scalar field, p B o and the cosmological constant A. Particularly, the study of this second 
stage involving essentially the two components, the BD scalar field $ and the cosmological 
constant A, can be thought of as our proposal for a k-essence model for dark energy. Then 
the difference between our dark energy model and those of other quintessence/k-essence 
theories lies in the fact that here in our model, we are interested in the way the presence 
of the BD scalar field, i.e., a k-essence, modifies the evolutionary behavior of the vacuum 
energy (A)-dominated epoch (which has an exponentially expanding nature in the standard 
general relativity context) while there, the quintessence/k-essence fields themselves (without 
the A-term) are expected to generate an accelerating expansion of some sort. 
The coupled Friedmann equation and the BD scalar field equation for the case at hand 
amount to setting p m = in eqs.(29) and (30), namely 




(44) 
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8tt 



,a/ (2a> + 3) 



(45) 



Again, these field equations should be supplemented by the remaining consistency equation 

A / <t> \ 

(46) 



PBD + H-) (PBD + P BD ) = Gq , $2 , • 



Again, we begin with the solution ansatz 



$(t) = 7r [l + X t]\ a(t) = a(0) [1 + X t] a . 



(47) 



Substituting these solution ansatz into the coupled field equations (44) and (45) yields in 
this time, 

(2o> + 1) 



a = 



7 = 2, 



(48) 



x= „ 2= 8^ , 2= (^ + 3 )( 6. + 5) 
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Thus the solutions behave as 



$(0 = i- [1 + X t] 2 , a(f) = a(0) [1 + X tf w+1)/2 ■ 

Having these exact solutions with us, again we consider the behavior of the BD scalar field 
(i.e., a k-essence) in this transition period. By plugging the exact solutions for the case at 
hand, eqs.(47) and (48) into eq.(ll), we are left with 



Pbd = -A 
Pbd = A 



4 (4a; + 3) 



a{t) 



-2/a 



(2cj + 3)(6c; + 5) \a(0) J 
8(3cu + 2) fa(t)Y 2/a 



(49) 



(2cj + 3)(6c; + 5) \a(0) / 



with a being given in eq.(48) and as a result, the equation of state of the BD scalar field 
perfect fluid in this transition period is given by 



Pi 



W B D = 



BD 



2(3u 



< 0. 



(50) 



Pbd (4a; + 3) 

In this transition period, however, the cosmic acceleration is supposed to set in. Thus one 
should demand a = (2a; + l)/2 > 1 or lu > 1/2 and this, in turn, leads to Pbd < 
and Pbd > as A > from eq.(48). Namely at this stage, the BD scalar field turns 
around and begins to behave as a positive pressure component making the total pressure 
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Ptot = Pa+Pbd = —A+Pbd /ess negative and hence slowing downihe accelerated expansion 
rate from the exponential law to the power law as we shall discuss in more detail below. 
Note, however, that since a(t) ~ t a and hence Pbd, Pbd ~ a~ 2 ^ a (t) ~ 1/t 2 , the positive 
contribution of Pbd to P to t = —A + Pbd would eventually become negligible as time goes on 
and the same argument holds for the total energy density p tot = Pa + Pbd = A + p B D namely, 
the negative contribution of pbd to ptot would become more and more negligible. This 
last observation, then, appears to indicate that the accelerated expansion would eventually 
"speed up" . Lastly, note that these solutions in eqs.(47) and (48) indeed satisfy the remaining 
consistency condition in eq.(46). It can be checked straightforwardly using the expressions 
for pbd and Pbd given in eq.(49) and the exact solutions in eqs.(47) and (48). Again, this 
consistency condition does not determine the BD ^-parameter to a particular value. 
These solutions in eqs.(47) and (48) have indeed been found some time ago in the so- 
called extended inflation model Q which attempts to resolve the graeeful exit problem and 
hence to save the original spirit of "old inflation" scenario. There the punchline was the 
realization that by replacing the general relativity content with the Brans-Dicke gravity, the 
accelerated expansion rate slows down from the exponential expansion to the power law one 
while the false-to-true vacuum phase transition rate remains the same and hence leading 
to the successful percolation and exit to the radiation-dominated era. Note again that in 
the standard Einstein gravity limit u —>■ oo in which the dynamics of the BD scalar field is 
washed out, i.e., $ — > 1/Go (as x "~ > 0), one recovers the familiar exponential expansion, 
a(t) = a(0)e*' in the Einstein gravity framework. Then the essential difference between 
the extended inflation then and the late-time cosmic acceleration now but with the same 
solutions would be ; (1) in the context of inflation at early universe, the model has to 
be constrained by such conditions as the enough expansion to resolve the difficulties with 
the hot big-bang model like the horizon, flatness and topological defects problems and the 
generation of density perturbation to some proper degree so as not to contradict with the 
observed anisotropy in CMBR (2) In the context of the late-time cosmic acceleration, 
on the other hand, the present dark energy model based on the BD theory should be tested 
carefully, say, by the current WMAP data [l( and the luminosity of type la supernova at 



high redshift 



12j both of which would provide a severe constraint on the current cosmic 



acceleration rate. Then eventually this would set a constraint on the value of the BD uj- 
parameter at this stage through the parameters a and \ in eq.(48). 
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3. Remarks 

(I) We go back to the case of MATTER-TO-SCALAR transition period studied above and 
examine the nature of solutions there. First, the condition of "cosmic expansion" amounts 
to demanding a = 2(u + l)/(3u + 4) > 0, which yields uo < —4/3 or w > —1. Secondly, 
there we demanded pbd > which led to w < —3/2 or —4/3 < uj < —6/5. Thus these two 
conditions, when combined, give u < —3/2. This, in turn, allows us to conclude that 



This inequality indicates that the expansion rate of the matter-dominated era in the presence 
of the BD scalar field, which is a k-essence, (or equivalently in the context of BD gravity 



of the k-essence field (or in the context of general relativity). Note also that from eqs.(35), 
(39) and (42), p m ~ a" 3 (t) ~ l/t 3a < p BD ~ a' 2 / a {t) ~ 1/t 2 as a > 2/3. Recall from 
section II that a(t) ~ t and thus pbd ~ a~ 2 (t) ~ 1/t 2 in the BD scalar field-dominated era. 
This observation indicates that as the evolution proceeds, the BD scalar field component 
takes over the energy density of the matter-dominated era and as a result, the BD scalar 
field-dominated era does indeed arrive as we assumed in the section II. And this supports 
our suggestion for the grand picture of universe evolution summarized in TABLE I. 
(II) Next, in the case of SCALAR-TO-ACCELERATION transition period, we learned that 
the accelerated expansion rate in the vacuum energy (A)-dominated era slows down to a 
power law by introducing the BD scalar field, which is again a k-essence (or equivalently in 
the context of BD theory) . 

Now these two observations (I) and (II) appear to imply that the effects of the presence 
of the BD scalar field (playing the role of an unique k-essence) on the late-time evolution 
of the universe can be summarized as follows ; the presence of the k-essence speeds up the 
expansion rate of the matter-dominated era (and hence reduces the deceleration) while slows 
down the expansion rate of the late-time cosmic acceleration phase (and hence again reduces 
the acceleration) ! This effect coming from the presence of the BD scalar field has been 
summarized in TABLE II. Thus to conclude, the BD scalar field, which is a k-essence field 
of our model, appears to interpolate smoothly between the decelerating matter-dominated 
era and the accelerating phase by reducing the deceleration of the former somewhat and 
then reducing the acceleration of the latter to some extent. And this role of "crossing 




(51) 



theory) namely, a{t) ~ t a , turns out to be greater than that, i.e., a(t) ~ t 2 l z in the absence 
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General Relativity (without BD scalar) 


Brans-Dicke Theory (with BD scalar) 


Matter-dominated era 


a(t) ~ t 2 / 3 


a(t) ~ t a , (a > 2/3) 


Accelerating phase 


a(t) = a(0)e** 


a(t) = a(0)[l + xt] (2w+1)/2 



TABLE II: Effect of the presence of the BD scalar field on the universe's late-time evolution. 



bridge" between the two late-time epochs is indeed consistent with the result of our earlier 
study of BD scalar field as an unique k-essence in section II in which we idealized the 
situation by ignoring all other types of matter, namely, the emergence of a zero acceleration 
epoch in between the two late-time eras. Thus we would like to point out that this smooth 
interpolation between the matter-dominated and the accelerating eras is the key role played 
by the BD scalar field as an unique k-essence. 

(III) It is interesting to note that the BD parameter u should "run" with the scale in order to 
serve as a successful model for dark matter - dark energy. Namely, according to the analysis 
performed in subsections 1 and 2 of section III above, it appears that the BD a>-parameter 
should behave as a running coupling parameter with growing behavior 

oo < -3/2 (say) in MATTER-TO-SCALAR transition period, 
< lo = —3/2 in (BD) scalar field-dominated era, 

> 1/2 in SCALAR-TO- ACCELERATION transition period 

as the scale factor a(t) increases (with time). As we mentioned earlier, the BD parameter 
uj in our model is a parameter of a k-essence (i.e., the BD scalar field) theory. Thus from 
a scalar field theory perspective, this "running" behavior of the BD w-parameter can be 
viewed as being natural in the sense of the renormalization group approach. What seems to 
be noteworthy is the fact that its value should grow with growing scale factor. We should, 
however, mention that this running behavior is not a prediction of our model. An explicit 
demonstration of this running behavior, although demanded, not only is beyond the scope of 
the present work but also seems unlikely to work out at the level of microscopic field theory 
as the present study involves the macroscopic equations of state in each epoch describing 
collective matter contents. 

(IV) Lastly, one may wonder if the presence of an additional epoch, i.e., the (BD) scalar 
field-dominated era could be inconsistent with the current estimate (and observation by 

n 

WMAP |1|) of the age of the universe based on the assumption that the universe is presently 
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matter-dominated and should have been so for most of its history. Under this assumption, 
using Hq = a/a = a/t, (a > 2/3) for the matter-domination and its observed value Hq ~ 
70km s' 1 Mpc' 1 with IMpc = 3.1 x 10 24 cm, the age of the universe is estimated to be 
to — 10 (Gyr). As is now well-known from, say, the WMAP data Q], the universe energy 
density consists of contributions of the same order of magnitude from dark matter and dark 
energy. And this indicates that the present stage of universe evolution is a sort of mixture of 
matter-dominated and accelerating eras. Furthermore, what we have realized in this work is 
that there could be an additional intermediate epoch of zero acceleration in between these 
two eras which is the (BD) scalar field-dominated epoch. Therefore, the conventional picture 
of universe history above now should be modified accordingly and the new ingredient realized 
in the present work, i.e., the (BD) scalar field-dominated era could also be inserted in the 
late time. And it means durations of all the three (allowing, of course, for their possible 
overlaps) would sum up to the age of the universe. To date, the cosmic age based on the 
best fit to a combination of WMAP, 2dFGRS (Two-Degree Field Galaxy Redshift Survey), 
Lya forest data, and the running index model is known to be to — 13.7 ± 0.2 (Gyr) |l|. 
In what follows, we attempt to work out this modification and to this end, we shall follow 
the formal way of estimating the age of the universe given by Peebles and Ratra Note 
first that even in the context of the BD theory, the Friedmann equation (in terms of the 
spatially-flat FRW metric) can be put in the form 



ay 8nG 



a J 3 
E(z) 



(prad + Pm + Pbd + Pa) = H E (z) where 



^r-ado(l + Zf + Q m0 {l + Zf + Q BD0 {1 + zf + tt A0 ] ^ , (52) 



"io = -TrjprPio, and (1 + z) - 



3# 2 r " ' a(t) 

with the subscript "0" indicating the present value of each quantity and z the redshift. That 
the Friedmann equation in BD theory can still be written as the first line in above equation 
(52) results from the fact that the terms u /6 — {a/a) say, in eqs.(36) and 

(44) can actually be rewritten in terms of (a /a) 2 using the solutions for the BD scalar field 
$(£) and the scale factor a(t). This then allows us to write the Friedmann equation in BD 
theory as in the first line in (52). Next, the age of the universe can be given by 

/"*« da _j f°° dz 
t0 = J T = H ° Lo (l + z)E(zY (53) 
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Although the formal evaluation of this integral is intractable, one may still wish to get a 
rough estimate of the cosmic age via the approximation 

rT , f f°° dz r z ™ dz 



\ka d (l + z)[Q rada (l + z)f /2 JzL (l + z)[fi m0 (l + ^) 3 ] 1/2 
r z BD dz r z i dz 1 

+ Un (l + z) in Bm (l + z)i] 1/2+ Jo (i + z) i^pj 

- H- 1 f-Q~ 1/2 1 I -Q~ 1/2 

1 2 [i + z j ,y 3 



+ o _1/2 



-rad) 

1 



:i + 4n) a + 4 



(l + z4)3/ 2 (1+^) 3 / 2 

+ fi A0 1/2 ln(l + z A )). 



The cosmological parameters in the standard LCDM model appearing in this expression for 
the cosmic age have been determined from the observations mainly by WMAP and they 
are O rad0 — 10~ 4 , fi m0 = VL vis0 + FIdmo — 0.05 + 0.20 = 0.25, f2 A o — 0.70 for the fractional 
present energy densities and (4 ad ~ 10 10 , ^ ~ 10 3 ), ~ 10 3 , 4 ~ 1), and 4 ~ 0.5 
(particularly note that the late-time cosmic acceleration had begun at the redshift of z ~ 0.5 
or 6 Gyr ago). The actual evaluation of the cosmic age in the context of our Brans-Dicke 
LCDM model based on eq.(54) above, however, is obscured as the data associated with the 
BD scalar field-dominated era such as VIbdo an d {z BD , z BD ) are unknown yet. In the most 
naive sense, however, we might wish to proceed with the rough estimate by assuming that 
since VL rad0 + Q m0 + &bdo + ^ao = 1> from the LCDM data above, perhaps flsDo < 0.05 and 
(z BD ~ 1, z BD ~ 0.5) as the BD scalar field-dominated era should be inserted in between 
the matter-dominated era and the current accelerating phase. Lastly, then, by plugging 
these data numbers in eq.(54) above, we end up with an estimate for the age of the universe 
in our Brans-Dicke LCDM model to be to — 1-63 x Hq 1 ~ 23 (Gyr). This is somewhat 
larger than the present estimate based on the best fit to observations. We hope, however, 
that more sensible choice of the data, fl B oo and (z BD , z BD ), if available at all, would bring 
it down to some extent. 

IV. Summary and discussion 

We now summarize our suggestion in the present work for adopting the BD theory of 
gravity as a "unified model" for dark matter - dark energy. First, it seems worthy of note 
that it has never been shown at the level of fundamental physics that the dark matter and 
the dark energy should be two distinct ingredients. Indeed in most of the recent studies on 
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the unified dark matter - dark energy, models employing a single component (whether it 
is the quintessence/k-essence or an exotic type of cosmic fluid) have been strongly favored. 
And needless to say, this is due to their attractive features of being able to provide an 
answer to the cosmic coincidence problem we mentioned earlier. In this regard, the way 
the BD theory serves as a unified dark matter - dark energy model is indeed quite different 
from that other models appeared in the recent literature do and thus is unique. Namely 
in our suggestion, first the BD scalar field can be thought of as being qualified to play 
a role of a k-essence since it possesses perhaps the simplest form of non-canonical kinetic 
term plus the interesting non-minimal coupling term to gravity. In section II, it has been 
demonstrated that the BD scalar field-dominated era is a zero acceleration epoch acting as a 
crossing bridge between the decelerating matter-dominated era and the accelerating phase. 
This result is remarkable enough but has been based upon the idealization/simplification 
that in the (BD) scalar field-dominated era, contributions from other types of matter may 
be dropped as they can be safely neglected. Therefore, in order to support and confirm this 
unique feature of the BD scalar field as a k-essence, more careful and realistic treatment is 
called for and this has been actually attempted in section III. That is, since the nature of the 
late-time universe evolution appears to be some mixture of matter-dominated, (BD) scalar 
field-dominated and accelerating phases, there we carried out a closer study of the effects of 
the BD scalar field on the evolutionary behavior of the matter-dominated and accelerating 
eras. And the result of the study in section III shows that the BD scalar field appears 
to interpolate smoothly between these two late-time stages by speeding up the expansion 
rate of the matter-dominated era somewhat while slowing down that of the accelerating 
phase to some degree. Thus with the newly found BD scalar field-dominated era in between 
these two, the late-time of the universe evolution appears to be mixed (allowing for possible 
overlaps) sequence of the three stages. Then our suggestion for employing the BD theory as 
a unified model for dark matter and dark energy can be stated as follows. In the MATTER- 
TO-SCALAR transition period, the dust-like matter and the BD scalar field together can 
be identified with some mixture of ordinary and dark matter and in the SCALAR-TO- 
ACCELERATION transition period, the BD scalar field and the cosmological constant A 
together may be thought of as constituting some exotic type of dark energy. In other words, 



unlike the other ambitious models in the recent literature 



13, 



14 Il5| attempting to build 



a unified dark matter - dark energy in terms of a single component, the philosophy behind 
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our model is to employ the BD scalar field as a catalyzer that essentially generates the dark 
nature from the ordinary matter and the cosmo logical constant. 

We now turn to the comparison of the feature our BD scalar field-dominated universe studied 
in section II with that of the curvature-dominated universe in the FRW universe model of 
Einstein gravity. As is well-known, the curvature-dominated universe is (effectively) an 
"empty space" solution in the spatially-open (k = —1) FRW model of Einstein gravity. As 
such, the universe expansion occurs essentially due to the negative curvature of the spatial 
section (k = —1). That is, all types of matter are absent (or assumed to be negligible) and 
the curvature term in the Friedmann equation behaves as if it is some kind of an energy 
density term, (a/a) 2 = —k/a 2 = 1/a 2 which gives a(t) ~ t. And the acceleration is zero 
simply because P = = p in eq.(21), i.e., (a/a) = — (4irG /3)(3P + p) = 0. 
By contrast, the BD scalar field-dominated universe can be thought of as the k-essence- 
dominated solution in the spatially-flat (k = 0) FRW model of a non-minimally coupled 
Einstein-scalar theory (or the pure BD theory of gravity). And the universe expansion 
occurs due to the k-essence (i.e., the BD scalar field) energy density pbd ~ 1/a 2 in the 
Friedmann equation, (a/a) 2 = (8iiGo/3)pbd = P 2 a 2 (0)/a 2 (where we used eqs.(5),(ll) and 
(15)). Particularly the zero acceleration can be attributed to the negative pressure Pbd — 
—Pbd/3 in eq.(24), i.e., (a/a) = — (8ttGo/9)(3Pbd + Pbd) = 0. To summarize, although the 
two universe models appear to exhibit the same expansion behavior a(t) ~ t, a — 0, the 
origin/nature of zero acceleration of the former comes from the negative spatial curvature 
(a geometrical nature) while that of the latter comes from the negative pressure of the BD 
scalar field or k-essence (a matter). Of course the essential difference between the two arises 
from the fact that the spatial section of the universe is open for the former model whereas 
it is flat for the latter model. 
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